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We present ground-state configurations and piiotoabsorp- 
tion spectra of Na-7+, Na-27+ and Na-41+. Botli tlie 
ionic structure and tiie pliotoabsorption spectra of medium- 
size sodium clusters beyond Na-20 have been calculated self- 
consistently with a nonspherical treatment of the valence elec- 
trons in density functional theory. We use a local pseudopo- 
tential that has been adjusted to experimental bulk properties 
and the atomic 3s level of sodium. Our studies have shown 
that both the ionic structure of the ground state and the po- 
sitions of the plasmon resonances depend sensitively on the 
pseudopotential used in the calculation, which stresses the 
importance of its consistent use in both steps. 
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Important progress has recently been made in the 
measurement of photoabsorjptijon spectra in small and 
medium-size sodium clusterscTu. This calls for a criti- 
cal re-evaluation of the theoretical description of the op- 
tical response of these systems. The earliest theoreti- 
cal invfajtigations of alkali clusters employed |t|he jellium 
modelu'EUmd were followed by other groupsBLi (for a re- 
view, seelia). The jellium model can qualitatively explain 
many of the experimentally observed featurestj. But a 
more realistic description of metal clusters is highly de- 
sirable, e.g. since the unrealistic sharp jellium edge gives 
rise to fundamental questions in the context of calculat- 
ing optical properties of metal clustersOo. 

The all-electron ab initio methods of quantum chem- 
istry treat the ionic degrees of freedota on the most so- 
phisticated level presently possiblellilll3. But their com- 
putation times grow so rapidly with the size that only 
small systems could be studied so far. Treating only the 
valence electrons explicitly and describing nucleus plus 
core electrons by a pseudopotential leads to a consid- 
erable simplification of the electronic system. The unre- 
stricted three-dimensional search foij-cluster ground-state 
configurations using Car-Parrinellot3't3 nn-Monte-Carlo 
methods with non-local pscudopotentialslij is, however, 
still a task of considerable complexity. 

Several approximate methods have been developed 
for the study of the ionic structure in larger clus- 
ters. T he || "s pherically averaged pseudopotential scheme" 
(SAPS)li30 optimizes the ionic positions in three dimen- 
sions while restricting the density of the valence electrons 
to spherical symmetry; the ion-valence electron interac- 
tion is described by the simplest pseudopotential possi- 
ble, the Ashcroft empty-core potential. However, where 
more accurate methods can be applied, these predict 



ground-state geometries which differ considerably from 
the SAPS results. The pseudopotential perturbation the- 
ory and similar approaches greatly improve on the SAPS 
deficiencies in the treatment of the valence electrons, but 
they either considjes only the volume-averaged effects of 
ionic i^mcturesBlla, or require their geometries as an 
inputCJ^EJ. Approaching the problem from the opposite 
direction, a sophisticated extension of the Hiickel modelcHi 
focuses on the prediction of ground-state configurations 
without explicitly taking the valence electrons into ac- 
count. Its results are in good agreement with ab initio 
calculations. However, the parameters of this model must 
be adjusted to ab initio calculations and the optical re- 
sponse of the electrons cannot be calculated. 

With the "cidindrically averaged pseudopotential 
scheme" (CAPS)a, a method has been developed that al- 
lows one to self-consistently calculate the ground states of 
clusters with several tens of atoms including ionic struc- 
ture without restricting the valence electrons to spherical 
symmetry. The ionic configuration is hereby optimized 
by the method of simulated annealing that seems to be 
the best method for coping with the strong isomerism 
found in larger clusters. The electmnic system is de- 
scribed in density functional theorytj; here we use the 
local density approxipaAtion (LDA) with the functional 
of Perdew and WangEEl. By an interlaced iteration, the 
set of equations 
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is simultaneously solved self-consistently. Here R denotes 
the set of all ionic positions, n the electronic density and 
E the energy functional 



E[n;R] = + E^c[n] + ^ J J d^r' d^r 

+ J n(r)Fei(r; R) d^r + ^ ^f.^i Tr-^ 

for a cluster of N ions with Z valence electrons each. 

The efficiency of the scheme results from two approx- 
imations that are made in the evaluation of the above 
energy functional. First, the interaction between valence 
electrons and ionic cores is described by a local pseudopo- 
tential: 



N 



yei(r;R) = VFps(|r-R, 



(1) 
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In the present work, we have developed a more physical 
pseud©Botential than that used in earlier applications of 
CAPSc^; it will be discussed below. 
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The second approximation is that while the ions are 
treated three-dimensionaUy, the electron density is re- 
stricted to cylindrical symmetry, i.e., in the solution of 
the electronic problem V"ps(|r — Rj|) is replaced by its 
cylindrical average 

1 

Vpsiz,p;z,,p,) = — Fp,(|r-R,|)d^. (2) 

This certainly is a simplification whose detailed conse- 
quences are hard to judge a priori. However, since the 
photoabsorption cross sections of singly charged sodium 
clusters show that the electronic density of mosl-daisters 
has an overall prolate, oblate or spherical shapeQ'EiEB, this 
approximation does not seem unreasonable for such sys- 
tems. A severe test for the quality of this approximation 
will be a comparison of its results to those of fully three 
dimensional methods. 

Before we present the results of our calculations, 
we discuss the local pseudopotential in some more 
detail. The most rigorous pseudopotentials in the 
sense of Philipps aad KleinmanE3 and modern ab ini- 
tio pseudopotentialsES are always non-local in the sense 
that each angular momentum component of the valence 
electron feels a different potential. However, the use of 
non-local pseudopotentials in the search for cluster con- 
figurations quickly exhausts computational resources be- 
cause of the multiple projections that have to be done at 
every step of the calculation, separately for each ion in 
the absence of any symmetry. Also, it has been shown 
that some ab initio pseudopotentials do nc^inecessarily 
lead to a good agreement with experimented. It there- 
fore makes sense to address the pseudopotential question 
from a more pragmatic point of view. 

Already early in the development of pseudopotential 
theory, it has been noted that by relaxing the Phillips- 
Kleinman condition, one can open up a new class of 
pseudopotentialstEl. They are called phcnomenological 
pseudopotentials or model potentials since they are con- 
structed by choosing some analytical function as a model 
potential and adjusting its parameters to experimentally 
known quantities, e.g., an atomic energy level or some 
bulk properties. Such model potentials can be nonlocal, 
or several-«artial-wave components may be chosen to be 
the sameE3. For metals with a simple electronic struc- 
ture like that of sodium, one can in this way construct 
pseudopotentials that are effectively local. Various lo- 
cal pseudopotentiak^ave, in fact, been successfully used 
(see, e.g., Refs.ESEilH). But the question of how a valid 
local pseudopotential should be constructed is nontriv- 
ial. In Ref.Ej this question was addressed in detail with 
an emphasis on solid-state properties. There, an "evanes- 
cent core potential" was proposed, and we have used it in 
some test cases. Most of our calculations, however, were 
done with a pseudopotential that we have constructed es- 
pecially for the use in finite sodium clusters, as explained 
below. 

Whereas in solid-state physics an important criterion 
for the practical usefulness of a pseudopotential is its 



fast convergence in reciprocal space, our aim is its ef- 
ficient handlipg in real space. We take up experience 
from CAPy2j and parameterize the pseudopotential on 
the basis of pseudodensities, related to (^) via Poisson's 
equation, which can be angle averaged analytically. This 
allows for a more efficient solution of the Coulomb prob- 
lem. Moreover, the short range of the pseudodensities 
ensures a fast repositioning of the ions. Our pseudopo- 
tential has the parameterization 

{^Qir^+ci r < ri 
-f r>r2. 

This corresponds to a pseudodensity with a two-step pro- 
file. Four of the seven parameters are fixed by requiring 
continuity of Vps{r) and its derivative. The remaining 
parameters gi , ri and r2 determine the physical proper- 
ties of the pseudopotential. Clusters contain from a few 
up to several thousand atoms, thus spanning the region 
from the atom to the bulk material. Our aim is therefore 
to develop a local pseudopotential that interpolates be- 
tween atomic and bulk properties. Thus we choose the 
parameters such that two quantities are reproduced cor- 
rectly: the atomic 3s energy level Cq on one hand, and the 
bulk Wigner-Seitz radius on the other hand. The lat- 
ter is determined by the minimum of the bulk energy per 
electron eb in second order perturbation theoryo. UsiM 
Ts = 3.93ao and the experimental value — —0.38 Ryu2l 
fixes Qi = — 0.503ao~'^ and r2 = 3.292ao and gives a con- 
straint on ri . The remaining freedom in choosing ri was 
exploited to fit the bulk compressibility B as closely 
possible to its experimental value Bcxp = 0.073 MbarES, 
yielding B = 0.0739 Mbar for n = 0.641ao. With these 
parameters we obtain Cb = —6.20 eV, close to the ex- 
perimental value —6.25 eV. The interstitial density, de- 
fined as the difference between the number of valence 
electrons in the Wigner-Seitz cell and in the muffin-tin 
sphereO, takes the value 0.35 with our pseudopotential. 
This agrees within 3% with the value given in Ref.E3. 
TheJaand-structure energy is 0.15 eV, in agreement with 
Ref.y. 

As a test-case study, we have calculated the ground 
state of Na-7+ with mw- pseudopotential, with the 
evanescent-corexifltentiaES, and with the empty-core like 
pseudopotentiaS. The £'5/i-geoniet£y of Na7-|- is well 
known from ab initio calculations£3ll3, and wc find the 
pentagonal bipyramid, shown in Fig. 1, with all three 
pseudopotentials. This demonstrates that CAPS can give 
realistic results even for very small systems, and that the 
cylindrical averaging is not too restrictive. The influ- 
ence of the pseudopotential can be seen in the bonding 
lengths, e.g. the distance between the two edges of the 
bipyramid: our pseudopotential and the evanescent core 
potential result in a distance of 6.03ao and 6.02ao, re- 
spectively, whereas the empty-core like pseudopotential 
leads to a shorter distance of 5.64aQ. All these values 
lie in the range found in fully three-dimensional calcu- 



2 



lations: Ref.Ej quotes a distance of 5.5ao, Ref.llll one of 
6.26ao. 
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FIG. 1. Right-hand side: CAPS-geometry of Na-7+. 
Left-hand side: percentages of the dipole mi sum rule ob- 
tained in the "localized" random-phase approximation for 
different pseudopotentials. Full line: present rtseudopo- 
tential; dashed line: evanescent-co*ei potentialEa; dotted 
line: empty-core like pseudopotential22l. A phenomenologi- 
cal Lorentzian line broadening is applied. 

Although these differences might appear smaU, they 
have noticeable consequences for the photoabsorption 
spectrum, as shown in Fig. 1. Here we plot the per- 
centages of the energy- weighted dipole sum rule (toi) 
obtained in rihe "localized" random-phase approxima- 
tion (LRPA)EfE3. For all pseudopotentials, two dominant 
peaks are seen that together carry about 90 % of the to- 
tal oscillator strength. Their relative heights reflect the 
oblate deformation of the electron cloud in the D^h con- 
figuration. (The remaining 10 % of the dipole strength 
are scattered in a region around-|4.5 eV.) The experi- 
mental photoabsorption spectrumH of Na7+ shows these 
two peaks at 2.35 eV and 3.25 eV when measured at 
a temperature of 650 K; at 35 K the lower transition 
splits into two close-lying peaks centered at 2.5 eV, while 
the higher transition is shifted to 3.3 eV. Our results do 
not reproduce the fine splitting of the lower peak, ft is 
knownE3 that the LRPA cannot resolve splittings that are 
due to interference of the plasmon resonance with specific 
particle- hole excitations, but it correctly reproduces the 
average peaks in the global strength distribution. (The 
convergence of the LRPA basis is within 3% for the peak 
positions.) The most important observation is, however, 
that all the spectra are blueshifted with respect to the ex- 
periment, and that the extent of the blueshift depends re- 
markably on the pseudopotential employed in the calcu- 
lation, ft ranges from 0.2 eV, found for the present pseu- 
dopotential, to nearly 0.5 eV for the empty-core type, 
with the evanescent-core results lying in between. That 
the empty-core like pseudopotential is the least accumte 
is not astonishing, since its parameters were fitteda to 
the bulk value of rg only in first-order perturbation the- 
ory. The discrepancy between the evanescent-core pseu- 
dopotential and ours could be the consequence of a bet- 
ter transferability of our pseudopotential; it might also 
reflect the fact that the evanescent-core pseudopotential 
is more difficult to handle numerically in our code. 



A blueshift of the dipole resonance in sodium clus- 
ters with respect to its experimental position has been 
found in many calculations, and thecR-lifis. f^lffi a long- 
standing discussion about its originotj'Ej'LjCj'Ej. Our re- 
sults show that the detailed form of the pseudopotential 
does strongly infiuence the resonance position, but that 
other effects must also contribute to the discrepancy with 
experiment. We think that one important effect is the 
finite temperature, present in most experiments, which 
results in an increased cluster volume, and thus to a de- 
crease of the plasnia frequency of a few percent. The 
experiments of Ref.cl show, indeed, that the average res- 
onance positions are shifted to lower energies when the 
temperature is increased. 




FIG. 2. Ionic structure and photoabsorption cross section 
a of Na-27-|- in arbitrary units against eV. Upper panel left: 
Experimental data£j. Lower panel left: LRPA photoabsorp- 
tion spectrum. Right-hand side: CAPS structure. 



This must be borne in mind when comparing the re- 
sults of calculations performed at T=0 to those experi- 
mental data of larger clusters, where the photoabsorption 
was measured only at finite temperature. On. the other 
hand, the LDA leads to^ell known errorscj that also 
affect the optic responseEa. 

Encouraged by the correct prediction of the structure 
of Na-7-t-, we have employed our scheme to calculate 
ground-state structures and photoabsorption spectra of 
sodium clusters in a size region where no self-consistent 
calculations with ionic structure have been made so far. 
In Figures |^ and |[ we present the ionic geometries and 
photoabsorption spectra of Na-27-|- and Na-41-|-. In both 
cases, we find good agreement with the experimental re- 
sults. The overall blueshift of the resonance peaks of 7 - 
9 % can again be accounted for by the effects mentioned 
above. The two pronounced peaks observed in Na-27-)- 
are a consequence of an overall prolate arrangement of 
the ions. Together they exhaust ^ 80% of the total os- 
cillator strength. The remaining strength is distributed 
around 4.3 eV, a region which was not scanned experi- 
mentally and thus omitted from the plot. 
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FIG. 3. Same as Fig. 1, but for Na-41+. Experimental 
data are from RefBI. 

In contrast, Na-41+ has a "magic" configuration of 40 
valence electrons, leading to a nearly spherical density 
and almost no splitting of the dipole resonance. The two 
strongest close-lying transitions exhaust 64% of the total 
oscillator strength. It is very interesting to note that, 
besides 20% of the strength lying outside the region of 
the experiment, another 16% of the strength are scat- 
tered around 3.4 eV. In the LRPA, these subpeaks are 
interpreted as internal compressional modes and modes 
of higher multipolarity that are coupled to the dipole os- 
cillation by the ionic structure. Their contributions are 
seen to be in good qualitative agreement with the high- 
energy shoulder observed in the experiment. 

In summary, our studies have shown that CAPS, com- 
bined with the localized random-phase approximation, 
provide a powerful tool for calculating ionic geometries 
and photoabsorption spectra of medium-sized sodium 
clusters where ab initio procedures are computationally 
not possible with a full search of ionic structures. With 
a local pseudopotential fitted to bulk and atomic prop- 
erties, we reproduce the global features of the dipole res- 
onances, apart from a slight blue shift that leaves room 
for temperature and non-LDA effects. We believe that 
CAPS gains additional importance since its results can 
serve as an input to three-dimensional calculations, thus 
bringing closer an explanation of the recent observation 
of the melting of clusters which depends on a detailed 
knowledge of their ionic structureO. 

This work has been partially supported by the 
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